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Abstract 

The energy dependence of the local V and CV violation in Au+Au and 
Cu+Cu collisions in a large energy range is estimated within a simple phe- 
nomenological model. It is expected that at LHC the Chiral Magnetic effect 
(CME) will be about 20 times weaker than at RHIC. In the lower energy 
range this effect should vanish sharply at energy somewhere above the top 
SPS one. To elucidate CME background effects a transport model including 
magnetic field evolution is put forward. 



1 Introduction 

As was argued in Refs. [HEllSlll] the topological effects in QCD with induced 
chiral asymmetry may be observed in heavy ion collisions directly in the 
presence of very intense external electromagnetic fields due to the "Chiral 
Magnetic Effect" (CME) as a manifestation of spontaneous violation of the 
CV symmetry. First experimental evidence for the CME identified via the 
observed charge separation effect with respect to the reaction plane has been 
presented by the STAR Collaboration [6]. In this paper we analyze the STAR 
data in a simple phenomenological way to estimate a possibility observing 
the CME in the larger energy range, from the LHC to FAIR/NICA energies. 
We also make a step towards a dynamical estimate of the CME background 
based on the nonequilibrium Hadron- String-Dynamics (HSD) microscopical 
transport approach [7] including magnetic field. 



2 Phenomenological estimates of the CME 

A characteristic scale of the process is given by the saturation momentum 
Qs DP) so the transverse momentum of particles pt ~ Qs- Then the total 
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transverse energy per unit rapidity at mid-rapidity deposited at the formation 
of hot matter is expressed through the overlapping surface of two colhding 
nuclei in the transverse plane S 

^ ~ e-V = e-Az.S = Q,-iQ',S)r^Q,-^^^^^ . (1) 
dy dy 

Here the energy density and longitudinal size Az ^ At ^ l/Q^ are taken in 
order of magnitude as follows e ~ and Az ~ l/Qs- 

For one-dimensional random walk in the topological number space the 
topological charge (winding number) generated during the time tb, when 
the magnetic field is present, may be estimated as 



= = V^TVV:^ ~ ^^^^^ (2) 

where is the sphaleron transition rate which in weak and strong cou- 
pling F5 ~ with different coefficients. The initial temperature Tq of the 
produced matter at time r ~ l/Qs is proportional to the saturation mo- 
mentum Qs, Tq = c Qs- At the last step of ([2]) the expansion time and 
the corresponding time dependence of the temperature are neglected. Since 
sizable sphaleron transitions occur only in the deconfined phase, the time tb 
in Eq. is really the smallest lifetime between the strong magnetic field tb 
one and the lifetime of deconfined matter r^: 

Tb = mm{fB,T,}. (3) 

The measured electric charge particle asymmetry is associated with the 
averaged correlator a by the following relation [8]: 

{C0s{i)a +i'(3- 2^i?,p)) = (cOs(V'c, + - '2lpc))/v2,c = Vi^aVl,(3 " 0^0/3 , (4) 

where \I/rp is the azimuthal angle of the reaction plane defined by the beam 
axis and the line joining the centers of colliding nuclei. Averaging in (jl]) 
is carried out over the whole event ensemble. The second equality in (jl]) 
corresponds to azimuthal measurements with respect to particle of type c 
extracted from three-body correlation analysis [8], vi and V2 are the directed 
and elliptic flow parameters, respectively. According to Ref. [Ij an average 
correlator a = ^Ja^^ is related to the topological charge. 
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where absorption and rescattering in dense matter responsible are neglected 
for the difference of magnitudes between the same and opposite charge cor- 
relations. In the last equality we assumed that Ql ~ s]^^ ~ dNh^dmns/dy [9]. 
Our susequent consideration is based on Eq. ([5]). 

Thus, the direct energy dependence of average correlator is comparatively 
weak. Results of dynamical heavy-ion calculations of the magnetic field at 
the central point of the transverse overlapping region of colliding nuclei and 
energy density of created particles are presented in Figs. [1] and [21 respectively. 
Here for a field estimate we follow Ref. [10] basing on the UrQMD model [11] 
and applying the electromagnetic Lienard-Wiechert potential with the retar- 
dation condition for the magnetic field. As is seen, at the impact parameter 
6 = 10 fm the maximal strength of the dominant magnetic field component 
By (being perpendicular to the reaction plane) is decreased in Au-|-Au colli- 
sions by the factor of about 10, when one proceeds from ^/sNN =200 GeV to 
Eiab =11 GeV, while for the created particle energy density e in the central 
box this factor is 250, i.e. noticeably higher. 
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Figure 1: The time evolution of the magnetic field strength eBy at the central 
region in Au+Au collisions with the impact parameter 6 = 10 fm for differ- 
ent bombarding energies. Calculations are carried out within the UrQMD 
model [11] (for a detail see [TO]). 

To use Eq. we need to identify the impact parameter, saturation 
momentum and multiplicity at a specific centrality. These can be found in 
Ref. [9] where the Glauber calculations were done. As a reference point we 
choose b =10 fm in our subsequent consideration. 

The measured value of {cos{ipa+4'i3—'^'^ rp)) for the same charge particles 
from Au-|-Au (a/JJvw = 200 GeV) collisions at the impact parameter b =10 
fm (40-50% centrality interval) is -(0.312 ± 0.027) ■ 10"^ [6j. Appropriate 
number for ^/snn = 62 GeV seems to be a little bit larger but for Cu-|-Cu 
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Figure 2: The time evolution of the energy density e of created particles 
in the Lorentz-contracted box with the 2 fm side at the central point of 
overlapping region. The impact parameter 6 = 10 fm. 



collisions the effect is definitely stronger [6]. Thus, ignoring any final state 
interactions with medium, assuming aa = = a and neglecting the directed 
flow Via = "Wife = we get from Eq. (jlj) a^^^ = 0.31 • 10^'^ for the maximal 
RHIC energy. Using numbers for the ^/snn =200 GeV reference case, from 
Eq. ([5]) we may quantify the CV violation effect by the correlator 

= Kau{V^)-^^^ -TB . (6) 

The normalization constant K^u can be tunned at the reference energy 
y^J/VTv =200 GeV from the inverse relation and experimental value a^xp at 
this energy for 6 =10 fm 

a2 . (200)1/8 

Kj^^ = -BP^ 1 — . (7) 

^" rB(200) ^ ' 

The lifetime tb may be defined as the time during which the magnetic field 
is above the critical value needed to support a fermion Landau level on the 
domain wall eBcru = 27r/S'rf, where Sd is the domain wall area. Since the size 
of the domain wall is not reliably known, it is hard to pin down the number, 
but it should be of the order of m^. Honestly, we have to treat it as a free 
parameter. 

Indeed the size of the topological defect (say, a sphaleron) in the region 
between Tc and 2Tc is very uncertain. At weak coupling, the size is deter- 
mined by the magnetic screening mass and it is ~ l/(asT). If one plugs 
as ~0.5 and T = 200 MeV, the size is of about 2 fm and then the threshold 
field is very small eBy ~ (agT)'^ ~ 0.2 m^. 
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Table 1: Estimated parameters for the CV violation effect in Au+Au colli- 
sions at centrality (40-50)% with the critical field eBcrit = 0.2 m^. 



^/SNN GeV 


1/16 


tb, fm/c 


Te, fm/c 


a' 


4.5 ■ 10^ 


2.86 


0.018 


>1 


0.016-10"^ 


200 


1.94 


0.24 


>2 


0.31-10-3 


130 


1.84 


0.33 


~2.3 


0.45-10-3 


62 


1.68 


0.62 


~2.2 


0.93-10-3 


17.9 


1.43 


1.41 


~2. 


2.48-10-3 


11. 


1.35 


1.66 


~ 1.9 


3.10-10-3 


4.7 


1.21 


0. 


0. 


0. 



On the other hand, we know that between Tc and 2Tc the magnetic screen- 
ing mass which determines the size of the sphaleron is not small as expected 
from the perturbative theory, agT, but from the lattice it is numerically large 
till about 5Tc. This would increase the threshold to 20 m^, however the re- 
lation between magnetic mass and the sphaleron size is valid only as long as 
the coupling is weak. All we can say it is perhaps in between (0.2—20) m^. 
Eventually lattice QCD calculations may clear this up. 

The upper bound on the magnetic strength eBcru = 20 results in Tb = 
even for the RHIC energy and therefore in this case the CME should not 
be observable at all in this energy range. The time evolution of the magnetic 
field and energy density, e, of newly created hadrons are presented in Figs. [1] 
and [21 The extracted values of tb defined by the constraints eBy > 0.2 
and r, (e > 1 GeV/fm3) are summed in Tabl. [H For the reference energy and 
the minimal magnetic field constraint we have Kau = 2.52 - 10-3. If lifetimes 
are known for all energies one can estimate the CV violation effect through 
the excitation function. 

From the first glimpse as follows from Tabl. [H in the case of eBcrit = 
0.2 the interaction time tb is defined solely by evolution of the magnetic 
field since tb < whereas ^ 2 fm independent of ^/SNN■ The expected 
CME for Au+Au at 6 = 10 fm (see the last column in Tabl. [1]) monotonously 
increases when y^J/vTv goes down but then sharply vanishes exhibiting a shal- 
low maximum in the range between near the top SPS and NICA energies. 
The position of CME maximum and its magnitude depend on the cut level 
which just defines tb- The decrease of the eBy bound till 0.02 shifts the 
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maximum toward lower energy ^/Jnn and enhances its magnitude. In an 
opposite limit when results are extrapolated to the LHC energy, the CME 
falls down by a factor of about 20 with respect to the RHIC energy. This 
result is quite understandable. The CME is mainly defined by the relaxation 
time of the magnetic field which is concentrated in the Lorentz-contracted 
nuclear region ~ 2R/'y. Therefore, the CME is inversely proportional to the 
colliding energy, ~ l/^SjVTv, and proceeding from the RHIC to LHC energy 
we roughly get the suppression factor about 4.5/0.2f» 22. 

There is one worrying point here. Proceeding from ^/SNN =200 to 62 GeV 
the predicted value of for 6=10 fm increases in three times though not more 
20% growth has been observed in these collisions in the recent experiment [6]. 
This essential disagreement cannot be removed by a simple variation of eBcru- 
One may try to explain this correlator overestimation at ^/sNN =62 GeV by 
an irrelevant choice of the energy dependence of multiplicity in Eq. ([6]). For 
the correlator ratio at these two energies we have 



where we use lifetime values from Tabl. [T] and experimental values for corre- 
lators [6], /3 = 1/8. As follows from Eq. ([8]), to explain the experiment the 
exponent should be negative, /3 < 0. Therefore, the fast growth of tb with 
the energy decrease cannot be compensated by uncertainty in the energy 
dependence of the correlator a. 

Uncertainty in the choice of the impact parameter does not help us to 
solve this issue. It turned out that one fails to fit this ratio by the variation 
of only eBcrit- Here we should remember that not only the strong magnetic 
field but also high density of soft equilibrium quark-gluon matter are needed. 
Equilibration requires some finite initial time tj^g which we associate with the 
moment when a maximum in the e{t) is achieved (see Figl2]). This makes 
f shorter and in combination with eBcrit variation, tb = tb — ti^s, allows 
us to satisfy the condition ([8]). Using the value of rB(200) obtained in this 
analysis one can recalculate the coefficient in Eq. K^u = 6.05 ■ 10~^, and 
therefore find the correlator a at any energy. In principle, similar analysis 
may be repeated for other impact parameters to consider the 6-dependence 
of the CME. As was shown in Refs. [3] the CME roughly is linear in h/R. 
Taking this as a hypothesis we evaluate the centrality dependence of the 
CME fitting this line to points 6 = 10 fm (or centrality (40 — 50)%) to be 
estimated in our model and 6 = where the CME is zero. The results are 
presented in Fig. [3] for Au+Au collisions at three energies. 




6 




Q. 















10 20 30 ^0 50 60 70 

% centrality 



Figure 3: Centrality dependence of the CME. Experimental points for 
Au+Au and Cu+Cu collisions are from [6]. The dotted line is our prediction 
for Au+Au collisions at the LHC energy. 

As it is seen the calculated lines quite reasonably reproduce the measured 
points of azimuthal asymmetry of charge separation for Au+Au collisions at 
^/sNN =200 and 62 GeV. The chosen value of eBcrit = 0.7 ml results in 
absence of the CME at the top SPS energy because the critical magnetic 
field practically coincides with the maximal field at this bombarding energy 
(see Fig. [T]). The CME at the LHC energy is expected to be less by a factor 
of about 20 as compared to that at the RHIC energy. Note that at the LHC 
energy we applied a simplified semi-analytical model [10] for magnetic field 
creation and assumed tj^g = 0. Thus, we consider this LHC estimate as an 
upper limit for the CME. 

Similar analysis can be repeated for Cu+Cu collisions basing on available 
RHIC measurements at two collision energies. Here one remark is in order. 
An enhancement of the CME in Cu+Cu collisions with respect to Au+Au 
ones was seen experimentally at the same centrality ^ but not at the same 
impact parameter. As follows from the Glauber calculations, the impact pa- 
rameter b=10 fm for gold reactions corresponds to centrality (40-50)% while 
the same centrality for copper collisions matches b=4.2 fm. The time distri- 
butions of the magnetic field and energy density for Cu+Cu collisions look 
very similar to that for Au+Au ones but lifetimes, both tb and r^, are shorter 
in the Cu+Cu case. For the extracted lifetimes and other characteristics at 
eBcrit = 0.2m^ {Kcu = 6.34 ■ 10^'^) we meet again the same problem: one 
should compensate a too strong energy dependence of the model correla- 
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tors by the proper definition of lifetimes. Defining tfie lifetime in the same 
manner as for Au+Au collisions the lifetime ratio tb{Q2)/tb{200) is turned 
out to be very close to experimental one at eBcrit =0.3 m^. In this case 
Kcu = 11-9 ■ 10~^. In the linear approximation with the reference point at 
b =4.2 fm, one may draw the centrality dependence of the CME for Cu+Cu 
collisions shown also in Fig|3] which is in a reasonable agreement with the 
experiment. Note that eBcrit =0.3 which is slightly above the maximal 
magnetic field at y^ijvw =62 GeV implies that the CME for Cu+Cu collisions 
will not be observable even at the top SPS energy. 

From dimensionality arguments the system-size dependence of the chiral 
magnetic effect (at the same all other conditions) would be expected to be 
defined by the surface S = S^ib) of an "almond" -like transverse area of over- 
lapping nuclei since both the high magnetic field and deconfined matter are 
needed for this effect. The magnetic field was evaluated in the the center of 
the overlapping region but as was shown in Ref. [10] the studied eBy com- 
ponent is quite homogeneous along x of this "almond" . Using for "almond" 
area a rough estimate as two overlapping discs of radius R = tqA^^^, namely 
S = SA{b) = n^R^ - (6/2)2(i? - 6/2), we have Scu{b = 4.2)/5au(& = 10) ^ 
1.65 which seems to be consistent with experimental ratio of the CME at 
y/SNN for these two points. However, this result was obtained for different 
eBcrit and non-zero initial time ti^^, and this success cannot be repeated for 
Cu+Cu (62 GeV) collisions. Therefore, the Cu enhancement effect is not 
only a geometric one. 

3 Towards a kinetic approach to the CME 
background 

The discussed CME signal may originate not only from the spontaneous 
local CP violation but also be simulated by other possible effects. In this 
respect it is important to consider the CME background. We shall do that 
considering a full evolution of nucleus-nucleus collisions in terms of the HSD 
transport model [7] but including formation of electromagnetic field as well 
as its evolution and impact on particle propagation. 

Generalized on-shell transport equations for strongly interacting particles 
in the presence of magnetic fields can be written as 

{§: + (Vp U) v.- - (v.- U + qvxB) V^^ } /(r,p, t) = Icoii{f, fi, ..Jn) 
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which are supplemented by the wave equation for the magnetic field whose 
solution in the semi-classical approximation for point-like moving charges is 
reduced to the retarded Lienard-Wiechert potential used above p[Q]. The 
term U ~ Re{Y7^^)/2pQ is the hadronic mean-field. 




Figure 4: Distribution of the magnetic field strength eBy in the y = plane 
at t =0.05 and 0.2 (in the middle) fm/c for Au+Au collisions at ^/snn =200 
and b =10.2 fm. The collision geometry is projected on x — z plane by 
points corresponding to a particular spectator position. Curves (and their 
projections) are levels of the constant eBy. 

One should note that the off-shell HSD transport approach is based not 
on the Boltzmann-like transport equation iQ but rather on the off-shell 
Kadanoff-Baym equations having similar general structure. The set of equa- 
tions was solved in a quasiparticle approximation by using the Monte-Carlo 
parallel ensemble method. To find the magnetic field a space grid was used. 
In a lattice point of this grid the retarded vector potential is evaluated. The 
magnetic field is calculated by its numerical differentiation. The field inside a 
cell is approximated by that at the nearest grid point. To avoid singularities 
and self-interaction effects, particles within a given cell are excluded from 
procedure of the field calculation. 

An evolution snapshot of the magnetic field By{x,y = 0,z,t) (in units 
of ml) formed in Au-|-Au (200 GeV) peripheral [b =10.2 fm) collisions are 
given in FigJUfor two time moments t =0.05 and 0.20 fm/c. The coUisional 
geometry is presented by a set of points every of which corresponds to a 
spectator nucleon. The whole field is not homogeneous exhibiting a wide 
maximum over the transverse size of overlapping (participant) matter and 
strong contraction in longitudinal direction. Opposite rotation of the mag- 
netic field along direction of two colliding nuclei results in corresponding two 
minima from outer sides of spectator matter remnants. At expansion these 
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remnants are moving away from each other. The position of a maximum 
in the magnetic field strongly correlates with that in the energy density of 
created particles. Large local values of By and e reached in these Au+Au 
collisions provide necessary conditions for observation of signals of a possible 
parity violation. 

4 Discussion and conclusions 

Summarizing one should note that for heavy- ion collisions at ^/SNN ^ H 
GeV the magnetic field and energy density of deconfined matter reach very 
high values which seem to be high enough for manifestation of the Chiral 
Magnetic Effect. However, these are only necessary conditions. To estimate 
a possible CME a particular model is needed. For the average correlator our 
qualitative prediction ~ ^n^n^ a rather small exponent but nevertheless 
it is too strong to describe observable energy behavior of the CME. This 
model energy dependence can be reconcile with experiment |6] by a detailed 
treatment of the lifetime taking into account both the time of being in a 
strong magnetic field and time evolution of the energy density in the QGP 
phase. For the chosen parameters we are able to describe data for Au+Au 
collisions on electric charge separation at two available energies. We predict 
that the effect will be much smaller at the LHC energy and will sharply 
disappear near the top energy of SPS. Coming experiments at the Large 
Hadron Collider and that the planned Beam Energy Scan program at RHIC 
[T2] are of great interest since they will allow one to test the CME scenario 
and to infer the critical magnetic field eBcru governing by the spontaneous 
local CV violation. 

The experimentally observed CME enhancement for Cu+Cu collisions 
is related with the selection of different impact parameters for the same 
centrality. However, it is not reduced to a purely geometrical effect. 

The problem of parity violation in strong interactions and the related 
CME are actively debated now. It is of great interest that the electric charge 
asymmetry with respect to the reaction plane may originate not only from 
the spontaneous local CP violation but also be simulated by other possible 
effects. First step in study of dynamical study of the CME background has 
been made in Sec|3l It is important that the developed kinetic approach in 
principle allows one to simulate the Chiral Magnetic effect itself. This work 
is in progress. 
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